On the strangeness content of the nucleon 
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We revisit the classical relation between the strangeness content of the nucleon, the pion-nucleon 
sigma term and the SU(3)f breaking of the baryon masses in the context of covariant chiral per- 
turbation theory. In particular, we consider the contributions of the decuplet resonances explicitly. 
We find that a value of the pion-nucleon sigma term of ~60 MeV is not at odds with, but favored 
" by the fulfillment of the Zweig rule. We compare these results with earlier ones and discuss the 

y—( _ convergence of the chiral series as well as the uncertainties of chiral approaches to the determination 

f^**; . of the sigma terms. 
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• The nucleon sigma terms, a^N and a s , are observables of fundamental importance which embody the internal scalar 

^-H ! structure of the proton and neutron. These are defined in the isospin limit (m u — m c i = rh) as 1 

S a„ N = -—(N\m(uu + dd)\N), 
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I. INTRODUCTION 



and contain information on the contribution of the quark-condensate to the masses of the baryons. Besides being an 
essential piece to the understanding of the mass of the ordinary matter, ovn and &s appear in the hadronic matrix 
elements of the (spin-independent) neutralino-nucleon elastic scattering cross section. An accurate determination of 
the sigma terms is then essential to constrain the parameter space of the underlying supersymmetric models from the 
experimental bounds in direct searches of weakly interacting dark matter particles [l| . 

It was first pointed out by Cheng that the contribution of the strange quark to the nucleon mass can be related 
^vq ■ with o^n and the SU (3) ^-breaking of the baryon masses in the octet [2:]. Namely, one can re-express the pion-nucleon 
sigma term as |3j 

C\l ■ a * N = i — ' \ 2 > 

. where y is the so-called "strangeness content" of the nucleon, 

• rH ■ 

X i 2{N\ss\N) = 2ma s = l _ 

■ (N\uu + dd\N) rrisG^N ct^jv' 

and 

TTl — 

a = ^ r (N\uu + dd-2ss\N). (4) 

Thus, (To can be understood as the value of the pion-nucleon sigma term in case that the strange-quark contribution 
to the nucleon wave function is null (Zweig rule 0)- These relations are very important in sigma-term physics as 
they imply strong constraints between o^n and y. In particular, requiring the Zweig rule to be approximately true 
(y ks 0) one obtains 



CttW ~ c . (5) 



1 We use the normalization of the nucleon spinors u(ff , s')u(p, s) = 2MmS s i s , with Mjy the nucleon mass. 
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Now, assuming linear SU(3)f breaking, one can write <j in terms of the physical masses of the baryon octet @ 
leading, in combination with Eq. ([5]), to 

m (M 5 + M s - 2M N ) 27 

o-ttW ; : MeV, (6) 

m s — to 1 — y 1 — y 

where we used m s /m — 26(4) [4]. The Eq. ^ is an estimate obtained at leading order in a S , J7(3)_F-breaking 
expansion of the baryon-octet masses. At first glance and given the success of the Gell-Mann-Okubo (GMO) relation 
3Ma + Me — 2Mn — 2Afs = 0, which is obtained at the same order in the expansion and is verified experimentally at 
the few-percent level, one expects this result to be a very good approximation. However, it is important to confront 
Eq. ([5]) with calculations of higher order terms in the SU (3) f -breaking expansion. 

This can be done in a model-independent fashion using baryon chiral perturbation theory (BxPT), which is the 
effective field theory of QCD at low energies @ and in the baryon sector [(| (for reviews see Refs. Q)- At leading 
order in the chiral expansion one recovers the linear 5[/(3).p-breaking corrections to the baryon masses and the GMO 
relation, whereas at next-to-leading-order (NLO), one obtains a contribution coming entirely from a chiral loop, so 
it is a pure prediction of B%PT (see e.g. the recent analysis in Ref. @). The NLO chiral correction to oo was first 
calculated by Gasser in Ref. 3] giving <tq = 35(5) MeV, although this calculation was done in the context of a chiral 
model of the meson cloud around the baryon which considered contributions from virtual octet baryons only (and 
not from decuplet resonances). Later, Borasoy et al. performed a calculation of the baryon masses and cjq in heavy- 
baryon (HB) 9] expansion and up to next-to-next-to-leading order (NNLO) 11] . In this work, the contributions of the 
decuplet resonances were not implemented explicitly but through resonance-saturation hypothesis they contributed 
to several of the many low-energy-constants (LECs) appearing at this order. All in all, they reported the value 
(To = 36(7) MeV, which was almost identical to the NLO result obtained by Gasser 15 years earlier. 

The agreement between these two different chiral calculations would suggest that the determination of gq in B^PT 
is under good theoretical control. Furthermore, its value is slightly below the canonical result a n N — 45 MeV [ldj . 
that was obtained by Gasser et al. from the irN elastic scattering data, 2 leading to a strangeness content in the 
nucleon in accordance with a rather small violation of the Zweig rule. This picture, however, has been challenged 
by recent studies based on modern ttN scattering data. Indeed, partial-wave analyses of the enlarged ttN database 



carried out by the George- Washington University group 14| , resulted in larger values of the pion- nucleon sigma term, 



&ttN = 64(8) MeV [15J. Besides that, a study of ttN elastic scattering in B%PT reveals that modern partial- wave 
analyses are, in general, more consistent with different scattering_phenomenology than the older ones and lead to a 
relatively large value of the sigma-term, cf. <t„n = 59(7) MeV [16]. These values of cr n N, together with the current 
one of <To(~ 35 MeV), imply a strangeness contribution to the mass of the nucleon of ^300 MeV. Although not 
impossible, such a scenario with a strong breaking of the Zweig rule is theoretically implausible, moreover after the 
experimental evidence pointing to a negligible strangeness contribution in other properties of the nucleon such as its 
electromagnetic structure [TtJ and spin (la |. 

In order to solve this puzzle the stress has been put, so far, on the value of the pion-nucleon sigma term. Inconsis- 
tencies and problems within the experimental irN database has been scrutinized (see e.g. Ref. [l9j]and the systematic 
errors in the relation between the isoscalar scattering amplitude and a^N carefully investigated [161 120( . Little atten- 
tion has been given to the value of ao though. There are two main lines along which the present determinations can 
be examined. On one hand, the result of Gasser Q can be put in a more evolved context than the one available at 
that time, whereas the one of Borasoy et al [ill ] might be afflicted by the poor convergence of the chiral series shown 
by HB in the SU(3)f theory p| [2l|. On the other, none of the two calculations include the explicit contributions of 
the decuplet resonances. These have to be introduced in SU (3) i?-BxPT due to the fact that the typical scale for the 
onset of the decuplet resonances (Mr — Mb) / A x s b ~ 0.3 is well below the chiral expansion parameter p = ttik / A x s b ■ 

A suitable approach to tackle both issues is a covariant formulation of B^PT with a consistent power-counting via 
the extended-on- mass-shell renormalization (EOMS) scheme (22|. The relativistic corrections included in this approach 
have been shown to tame t he p oorly convergent series of the HB expansion in baryonic observables as important as 
the magnetic moments [2l|, |23| or masses [1, [24|. Moreover, once a prescription is taken to treat the problem of the 
interacting Rarita-Schwinger fields [25[ , this scheme is straightforwardly applicable to include the contributions of the 
decuplet resonances [23| • In this work we calculate cto up to NLO using covariant BxPT renormalized in the EOMS 
prescription and including explicitly the effects of the decuplet. We will compare the results with those obtained in 
the HB expansion and we estimate systematic higher-order effects through a partial calculation of NNLO pieces. All 
together, we will find that the effect of the decuplet is sizable and that it pushes the value of ctq upwards making the 



2 For a detailed exposition of the dispersive methods for obtaining ct^jv from the analytic continuation of the ttN scattering amplitude to 
the Cheng-Dashen point see Refs. [l(tll2l [T3l . 
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(a) (b) (c) 

FIG. 1: Feynman diagrams contributing to the nucleon mass up to 0(p 3 ) in BxPT. The internal solid lines correspond, in 
general, to any octet-baryon, double lines to decuplet-baryons and dashed lines to mesons. The black dots indicate l s '-order 
couplings while crosses are insertions of 0(p 2 ) operators given by the LECs bo, bn and 6f responsible for the leading SU(3)f 
breaking of the baryon-octet masses. 



modern experimental determinations of a^N ~ 60 MeV consistent with the Zweig rule. A first indication that the 
decuplet contributions could help to solve the strangeness puzzle was given by the HB calculations in Ref. [26] and, 
indirectly, in Ref. [27j where very large and negative values of o s were obtained when demanding o^n = 45 MeV, 
indicating a larger ctq. 



II. CALCULATION 

The expressions for the sigma terms can be obtained either from the explicit calculation of the scalar form factor 
of the nucleon at q 2 = or applying the Hellmann-Feynman theorem to the chiral expansion of its mass, 

A dM N ml ( 1 1 1 \ 4 , 

**» = m ft*" = -f + 2^8^ + 3^0^) Mn + °^ 

m dM N 2 ml f 1 t 2 \ M + 4 

s s dm s K 2 \ 2mA' 3m^ 0m,, / 

We will follow the latter strategy since the explicit expressions for the baryon masses in the different schemes treated 
in this paper can be directly obtained using the Appendix of Ref. j8j. Thus, the chiral expansion of the sigma terms 
up to NLO from Eq. (J7J is written as, 

2 

m 

cr^N = -4(2& + b D + b F )-^- + 
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m± 



a s = -4(26 + b D - bp) [m z K 

jtJ^ E (4 B j4 B) (^)+4 T i^ T) (^))+o(/)- (8) 

^' <f>=n,K,ri 

The first line in these formulas corresponds to the LO contribution given at tree-level by the same 0(p 2 ) LECs that 
appear in the chiral expansion of the baryon masses. While &o provides a SU(3) ^-singlet contribution, the LECs bo 
and bp induce an octet breaking (tree-level in diagram (a) in Fig. [T]) which leads to a description of the baryon-octet 
mass splittings that is equivalent to the GMO relation [8]. The second lines enclose the NLO corrections coming, at 
0(p 3 ), from the loop topologies shown in Fig.Q] (b) and (c). Thus, the effect of virtual octet (B) and decuplet (T) 

(X) 

baryons is explicitly accounted for. Their contributions are weighted by the coefficients £ N L which are combinations 
of 517(3) Clebsch-Gordans and the meson-baryon couplings D, F (octet contributions) and C (decuplet contributions). 

(X) 

The loop functions Ea depend, exclusively, on the mass of the virtual pseudoscalar meson and on the ones of the 
octet and decuplet in the chiral limit, Mb and Mj respectively. Strictly speaking, the 5'C/(3)_F-breaking of the baryon 
masses in these loops, which are represented by the crosses in Fig. [1] (b) and (c), are contributions that start at 
NNLO or 0(p 4 ). 

For the baryon masses we use the results obtained in Ref. [8j in Lorentz covariant B%PT up to 0(p 3 ). The 
chiral loops contain divergences and analytic pieces breaking the power-counting formula [6j that are removed in 
dimensional regularization by the proper redefinition of the bare LECs (EOMS scheme [22]). The contributions of 
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the decuplet baryons are included taking the octet and decuplet masses in the chiral limit of approximately of the 
same order. Namely, the octet and decuplet contributions are considered on the same footing for power-counting 
purposes and no specific expansion in S = (Mt — Mb) is performed. The HB formulas (2(| can be always recovered 
from the renormalized covariant results by taking the non-relativistic expansion Mb ~ Mt ~ A x g#. In particular 
for the decuplet, the HB results within the small-scale-expansion (SSE) |28[ (which furthermore considers 6 ~ p) are 
recovered 1271 1 . 



TABLE I: Values of the 0(p 2 ) LECs bp, and bp determined from the baryon octet mass splittings in the different BxPT 
approaches considered in this paper. 





Tree level 0(p 2 ) 


Octet 0(p 3 ) 
HB Covariant 


Octet+Decuplet 0(p 3 ) 
HB-SSE Covariant 


b D [GeV" 1 ] 
b F [GeV" 1 ] 


0.060(4) 
-0.213(2) 


0.061(4) 0.061(4) 
-0.502(2) -0.420(2) 


0.315(4) 0.161(4) 
-0.704(2) -0.502(2) 



For the numerical values of the couplings, we use D = 0.80 and F = 0.46 [2^. The decuplet coupling C can be 
fixed from the A(1232) — > irN decay rate, giving C = 1.0 [23| . However, there is some evidence from lattice QCD 
that this coupling is somewhat smaller (24[. Indeed, an 5[/(3)F-average among the different decuplet-to-octet pionic 
decay channels gives C = 0.85 ± 0.15, that is the value we use. 3 As mentioned above, the 0(p 2 ) LECs bu and &f 
are determined using the experimental baryon-octet mass splittings. Their values for the different BxPT schemes 
analyzed in this paper can be found in Table U For the meson decay constant we also take the S£/(3).F-average 
F<j, = 1.17/tt with f„ — 92.4 MeV. For the masses of the pseudoscalar mesons we use m w = m w ± = 139 MeV, 
m,K = rriK± = 494 MeV, while for the baryon masses in the loops we use the chiral-limit baryon masses obtained at 
LO, M B 1] = 1.151 GeV and M T 1] = 1.382 GeV. The mass of the 77 meson is fixed with the Gell-Mann-Okubo mass 
relation, 3m 2 — Am 2 K — which is accurate enough up to the order we work. 

Finally, we will restrain our analysis to 0(p 3 ) despite of the fact the extension of formulas to 0(p 4 ) accuracy is 
straightforward and have been reported in the literature [TTI. l30l - [32l | . At the latter order, 15 new LECs contribute to 
the baryon masses and sigma-terms. Ei ght of them correspond to 0(p 2 ) operators which appear through diagrams 
with the topology of a tadpole (see Ref. [30] for details). These also contribute to the chiral expansion of the meson- 
baryon scattering amplitudes, although their LECs have not been determined yet from the associated experimental 
data or LQCD results. The other loop diagrams appearing at this order are the ones at 0(p 3 ) but with the SU(3)f 
breaking of the baryon masses in the loop taken into account by insertions of the 0(p 2 ) LECs bo, bu and bp (crosses in 
the diagrams (b) and (c) of Fig. [T]). The remaining 7 LECs correspond to 0(p 4 ) operators and they renormalize the 
loop divergences appearing at this order. Therefore, a quantitative analysis of the sigma terms at NNLO without any 
further assumption on the values of the LECs (such as Large JV C constraints [3l[ or resonance saturation hypothesis 
estimates [ill ]) is affected, at present, by a large uncertainty. On the other hand, a promising source of theoretical 
information on the values of the LECs is becoming available through lattice QCD calculations. An application in this 
direction and performed in the EOMS scheme at 0(p 4 ) can be found in 32]. 

Nevertheless, the analysis of part of the 0(p 4 ) corrections can be useful to asses the convergence of the chiral series 
and to give a credible estimate on the systematic error to the 0(p 3 ) results on the sigma terms due to the truncation 
of their chiral expansions. Indeed, we have calculated explicitly the respective corrections arising from the SU(3)p 
breaking of the baryon masses in the loops (b ) and (c) of Fig. [TJ The divergences have been renormalized in the 
EOMS scheme and the uncertainty on the unknown values of the 0(p A ) LECs has been explored by varying the 
renormalization scale in the interval 0.7 < /1 < 1.3 GeV. The maximal contribution obtained for these corrections in 
both, the octet and decuplet diagrams, is quoted as our theoretical uncertainty. That is, 

Ac™ ~ 20 MeV, Act™ ~ 140 MeV, 

Act^° ms ~ 6 MeV, Aaf OMS ~ 60 MeV. (9) 

This explicit calculation of higher-order pieces already confirms the expectation that the convergence in the covariant 
approach is substantially better than the one obtained in the HB approach [§|, [2l| . 



3 Note that the value for C of the present work is different from the one often used in HB calculations @. In these papers, a convention 
for the "vielbein" that is related to ours by a factor of 2 is employed. Moreover, the value we use in this paper is different to the one 
used in [j|, explaining the slightly different decuplet results obtained here and there. 
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III. RESULTS 

From the discussion above and the Eqs. (|SJ), it is clear that the only unknown parameter in the chiral expansion of 
the sigma terms up to NLO is the LEC bo . It accompanies a singlet structure that cannot be disentangled from the 
baryon mass term by looking at the experimental baryon masses alone ||. The strategy followed by us in relation 
with Eq. ((2|) is that one sigma-term can be predicted if the other is taken as input. Namely, ctq is the pion-nucleon 
sigma term once bo is fixed by the requirement a s — 0. And inversely, a s and the strangeness content of the nucleon 
can be predicted using the determination of ovjv obtained from the irN scattering data. In Table HI] we show the 
results for ao and 60 given by the exact fulfillment of the Zweig rule and for the different B^PT schemes treated in 
this work. The errors in bo are obtained assuming an uncertainty in the Zweig rule of the order of the systematic 
error estimate of Eqs. ((9|) for a s - This error propagates into ao and is added in quadrature with the systematic one 
in ovjv- 

TABLE II: Values of the 0(p 2 ) LEC 60 given by the fulfillment of the Zweig rule and the corresponding values of ao for the 
different B^PT approaches considered in this paper. 





Tree level 0(p 2 ) 


Octet 0{p 3 ) 
HB Covariant 


Octet+Decuplet 0(p 3 ) 
HB-SSE Covariant 


(Zweig [Gc y-l] 


-0.274 


-0.90(15) -0.70(5) 


-1.52(15) -0.95(5) 


ao [MeV] 


27 


58(23) 46(8) 


89(23) 58(8) 



As we can see in the Table [TT1 the corrections to the LO result on a studied are large. This occurs despite that the 
discrepancy of the GMO equation is correctly predicted in any of these schemes and, in fact, the description of the 
experimental octet mass splittings improves at 0{p 3 ) Q. In other words, the empirical success of the GMO equation 
does not preclude Eq. ^ of receiving sizable SU(3)p breaking corrections. Another important outcome concerns the 
comparison between the results obtained in HB and in covariant B%PT. As already anticipated by the calculation 
of the 0(p A ) pieces in Eqs. ©, the SU(3)f HB expansion has severe problems of convergence in the description of 
the sigma terms at 0(p 3 ). The huge central value and errors of ao for the HB-SSE expansion has to be regarded 
as a clear manifestation of these problems. Another indication in this sense is the large change in the values for ao 
between HB and HB-SSE in Table HTl a change 3 times larger than for the covariant calculation. Such large variation 
is at odds with the idea of resonance saturation of the LECs at 0(p 2 ), indicating that higher order contributions 
from resonance exchanges beyond 0(p 3 ) are still sizable. On the contrary, for the covariant calculation offered the 
difference between the calculations including/excluding the explicit decuplet of baryon resonances is only of around 
10 MeV, much smaller than the difference between the LO and NLO results. This indicates a stabilization of the 
final outcome at the 0(p 3 ) value for the covariant case. These conclusions are consistent with those derived from the 
analyses of other observables [3, Hi], [23| that also indicate similar problems of convergence for the HB studies in the 
SU(3)f sector. Hence, in the rest of the paper we focus on the results obtained in covariant B%PT. 

The result with only the octet contributions ought to be compared with the classical result of Gasser Notice 
that at the time this calculation was published xPT was m an early stage and not fully developed. Interestingly 
enough, this author noticed back then that using the leading non-analytic contribution (or LNAC, the old-fashioned 
label for the leading HB loop terms) to the baryon masses led to a very poor convergence of the respective chiral 
series. In order to overcome this problem, he proposed a relativistic and chiral model of the meson-baryon interaction 
and computed the non-analytic chiral corrections by the explicit calculation of the corresponding loop diagram, Fig. [T] 
(b). His approach was close, in spirit, to the one used in this paper although there are some technical differences. 
The most important one is that a cut-off was used for the regularization of the loop divergences. His final result, 
(To = 35(5) MeV, is obtained for a given value of the cut-off and it is not clear how much of the quoted uncertainty 
originates from reasonable variations of this parameter of the model. With these considerations, we conclude that the 
result obtained in our paper for the octet contribution is in reasonable agreement with the one of Gasser. 

TABLE III: Value of the LEC bo and of the observables related to the strangeness of the nucleon, y and er s , obtained in covariant 
BxPT including decuplet contributions and using the experimental determinations of ovat as input. 





C Pt ' \G0V- 1 } 


y cr s [MeV] 


(7-nN = 45(7) MeV 


-0.79(9) 


-0.28(13)(10) -150(80)(60) 


a^N = 59(7) MeV 


-0.97(9) 


0.02(13)(10) 16(80)(60) 



The new contributions given by the decuplet baryons are not negligible producing a ^10 MeV rise on ao. This result 
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is very important for sigma-term physics and it has a strong impact on statements about the strangeness content 
of the nucleon based on the value of o^n- In order to illustrate this point let us inverse the logic used before and 
determine the values for y and a s obtained from the two distinct experimental determinations ct^n — 45 MeV [l"ol | 
and ct„n ~ 59(7) MeV [161 . As we can see in Table IIII1 the determination of the strangeness in the nucleon depends 
strongly on the value of the pion- nucleon sigma term due to a relative factor ~ m s /(2rh) ~ 13 between the two 
observables. Hence, a variation of 10 MeV in ovjv causes changes in a s spanning more than 100 MeV. This factor 
also amplifies the uncertainty on the latter observable that propagates from the relatively small error of the former. 
One has to stress that this is a general feature relative to the different quark mass entering in the definitions of the 
sigma terms, Eqs. ([1]), which makes it very difficult to give predictions of a s with the ballpark accuracy of few tens of 
MeV. Two errors are quoted for y and a s in the last two columns of Table IIIII The first stems from the propagation 
of the uncertainty in <t„n and the last from the estimated 0(p 4 ) uncertainty in our calculation, Eq. ([§]). 

Nevertheless, the most important result in Table IIIII is that a pion-nucleon sigma term of ~ 60 MeV is perfectly 
consistent with an approximate fulfillment of the Zweig rule and, hence, to a small strangeness content in the nucleon. 
On the other hand, the canonical value of ovjv gives a prediction for the strangeness that is quite large and negative. 
It is then clear that the present calculation of do based on covariant B^PT in EOMS renormalization scheme with 
the decuplet as explicit degrees of freedom clearly favors relatively large values of o^n with small strangeness content 
in the nucleon. 



IV. CONCLUSIONS 



In summary, we have revisited an old empirical relation between the strangeness content of the nucleon and the 
pion-nucleon sigma term in the context of covariant B^PT. Earlier estimates of <To made at different levels of accuracy 
in xPT agreed on that a small violation of the Zweig rule in the nucleon requires a value of a^N close to ^35 MeV. A 



long-standing puzzle [33j has arisen from sustained experimental evidence pointing to a value of this quantity close to 
60 MeV, reinforced by the values obtained from modern ttN databases. We have shown that the previous calculations 
of (To are afflicted by important systematic effects, in particular those given by the omission of the decuplet resonances. 
Once these are incorporated, the relatively large value of a n N is consistent with and favored by the fulfillment of the 
Zweig rule. 
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